Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 2 February 2008 



(MN M?eX style file vl.4) 



Reionization History from Coupled CMB/21cm Line Data 



R. Salvaterra^, B. Ciardi^, A. Ferrara^, C. Baccigalupi 

^ SISS A/International School for Advanced Studies, Via Beirut J^, 34014 Trieste, Italy 
^ Max-Planck-Institut fiir Astrophysik, Garching, Germany 



2 February 2008 



in 
o 
o 

(N 
X) 

(N 



> 

o 
in 
o 

6 



ABSTRACT 

We study CMB secondary anisotropies produced by inhomogeneous reionization by 
means of cosmological simulations coupled with the radiative transfer code CRASH. The 
reionization history is consistent with the WMAP Thomson optical depth determina- 
tion. We find that the signal arising from this process dominates over the primary 
CMB component for I 4000 and reaches a maximum amplitude of l{l + l)C;/27r ~ 
1.6 X 10""'^'^ on arcmin scale, i.e. I as large as several thousands. We then cross-correlate 
secondary CMB anisotropy maps with neutral hydrogen 21cm line emission fluctua- 
tions obtained from the same simulations. The two signals are highly anti-correlated 
on angular scales corresponding to the typical size of Hii regions (including overlap- 
ping) at the 21cm map redshift. We show how the CMB/21cm cross-correlation can 
be used to: (a) study the nature of the reionization sources, (b) reconstruct the cosmic 
reionization history, (c) infer the mean cosmic ionization level at any redshift. We 
discuss the feasibility of the proposed experiment with forthcoming facilities. 

Key words: galaxies: formation - intergalactic medium - cosmology: theory 



1 INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP*) has 
provided strong evidence for an optical depth to Thomson 
scattering of Te 0.17±0.04 (the uncertainty quoted for this 
number depends on the analysis technique employed), based 
on the measured correlation between Cosmic Microwave 
Background (CMB) temperature and polarization on large 
angular scales (e.g. Kogut et al. 2003). If the reionization 
process is described as instantaneous and homogeneous, this 
corresponds to a reionization redshift Zion ~ 17 ± 5. More 
probably, reionization went through a highly inhomogeneous 
phase (e.g. Ciardi et al. 2000; Gnedin 2000; Miralda-Escude, 
Haehnelt & Rees 2000; Ciardi, Stoehr & White 2003, CSW; 
Sokasian et al. 2003; Ricotti & Ostriker 2004), which ended 
only when the individual Hii regions overlapped completely. 
In this case, the reionization process should have left an im- 
print on the CMB. In fact, the modulation of the ionization 
fraction, playing a similar role to the density modulation 
from the non-linear Vishniac effect, leads to anisotropies at 
sub-degree scales (e.g. Bruscoli et al. 2000; Benson et al. 
2001; Gnedin & Jaffe 2001; Santos et al. 2003). In addition 
to temperature anisotropies, Thomson scattering introduces 
a polarization signal in the CMB spectrum. The detection 
of anisotropies in the temperature/polarization power spec- 
trum is an invaluable tool to discriminate between different 
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sources of ionizing photons and reionization histories (e.g. 
Bruscoli, Ferrara & Scannapieco 2002; Holder et al. 2003; 
Naselsky & Chiang 2004). 

An alternative way to probe the end of the cosmic 'dark 
ages' is through 21cm tomography. From the pioneering 
work of Field (1959), it has been suggested that the neutral 
hydrogen in the Intergalactic Medium (IGM) and in gravi- 
tationally collapsed systems may be detectable in emission 
or absorption against the CMB at the frequency correspond- 
ing to the redshifted 21cm line associated with the spin-ffip 
transition of the hyperfine levels of neutral hydrogen. The 
inhomogeneities in the density field, ionized hydrogen and 
spin temperature produce signatures both in the angular 
and in the redshift space. Different signatures have been in- 
vestigated, ranging from the 21cm line emission induced by 
the 'cosmic web' (Madau, Meiksin & Rees 1997; Tozzi et al. 
2000), the neutral hydrogen surviving reionization (e.g. Cia- 
rdi & Madau 2003; Furlanetto, Sokasian & Hernquist 2004; 
Furlanetto, Zaldarriaga & Hernquist 2004) or the minihalos 
with virial temperatures below 10"* K (e.g. Iliev et al. 2002), 
to the 21cm lines generated in absorption against very high- 
redshift radio sources by the neutral IGM (Carilli, Gnedin & 
Owen 2002) and by intervening minihalos and protogalactic 
disks (Furlanetto & Loeb 2002). 

In this paper, we compute the CMB temperature 
anisotropies due to an inhomogeneous reionization history 
obtained from radiative transfer simulations consistent with 
WMAP observations (Ciardi, Ferrara & White 2003, here- 
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after CFW). Moreover, we cross-correlate them with the ex- 
pected 21cm emission maps obtained by Ciardi & Madau 
(2003, hereafter CM) for the same simulations, and discuss 
how the cross-correlation can be used to reconstruct the 
reionization history and to constrain the nature of ionizing 
sources. Our work is similar in spirit to the recently pub- 
lished study by Cooray (2004), although that work is based 
on a simplified analytical description of the reionization pro- 
cess. This might be the reason for which our conclusions 
differ from those obtained by Cooray (see Section 6). 

The paper is organized £is follows. In Section 2 we 
present the numerical simulations of IGM reionization by 
CSW and CFW, and in Section 4 the results of CM on the 
21cm emission from such patchy reionization histories arc 
briefly described. In Section 3 we construct and study the 
maps and the angular power spectra for secondary CMB 
temperature anisotropics due to the above reionization pro- 
cess, whereas the cross-correlation between CMB and 21cm 
maps is presented in Section 5. Finally, in Section 6 we sum- 
marize and discuss the results. 

Throughout the paper we adopt the ACDM "concor- 
dance" model with n^=0.3, nA=0.7, h =0.7, ^6=0.04, n=l 
and (78=0.9, within the WMAP experimental error bars 
(Spergel et al. 2003). 



2 NUMERICAL SIMULATIONS OF IGM 
REIONIZATION 

In this Section we briefly describe the numerical simulations 
of IGM reionization adopted to model the 21cm line emission 
from neutral IGM and the CMB temperature anisotropics, 
and refer to CSW and CWF for further details. 

A cosmological volume of comoving side 479/i.~^ Mpc 
has been simulated (Yoshida, Sheth & Diaferio 2001) with 
the N-body code GADGET (Springcl, Yoshida & White 2001). 
An approximately spherical region with a diameter of about 
50h~^ Mpc has been subsequently "re-simulated" at a higher 
resolution (Stoehr 2003) with the technique described in 
Springel et al. (2001, hereafter SWTK). A friends-of-friends 
algorithm was employed to determine the location and mass 
of dark matter halos. Gravitationally bound substructures 
have been identified within the halos with the algorithm 
SUBFIND (SWTK) and have been used to build the merg- 
ing tree for halos and subhalos following the prescription of 
SWTK. A particle mass of Mp = 1.7 x W^h~^ Mq allows 
to resolve halos as small as M ~ 10^ M©. The galaxy pop- 
ulation has been modeled with the semi-analytic technique 
described in Kauffmann et al. (1999) and implemented as in 
SWTK. For each of the simulation output we compile a cat- 
alogue of galaxies containing for each galaxy, among other 
quantities, its position, mass and star formation rate. 

A cube of comoving side L = 20h~^ Mpc has been cut 
from the high resolution spherical subrcgion to model the 
details of the reionization process, using the radiative trans- 
fer code CRASH (Ciardi et al. 2001; MascUi, Ferrara & Ciardi 
2003) to follow the propagation into the IGM of the ionizing 
photons emitted by the simulated galaxy population. Several 
sets of radiative transfer simulations have been run in CSW 
and CFW, with different choices for the gala^xy emission 
properties. The ones used here are labeled S5 ('late' reion- 
ization case) and L20 ('early' reionization case), and adopt 



an emission spectrum typical of Pop III stars, a Salpeter Ini- 
tial Mass Function (IMF) and an escape fraction of ionizing 
photons /esc = 5% (S5) and a Larson IMF with fesc = 20% 
(L20). For details and discussion on the choice of parameters 
we refer to CSW and CFW. The S5 and L20 simulations give 
a reionization redshift of Zion ~ 8 and « 13, respectively. In 
addition, they provide the redshift evolution and the spatial 
distribution of ionized/neutral IGM and have been used to 
model both the Hi 21cm line emission (see Sec. 4) and the 
CMB temperature anisotropics (see Sec. 3). 



3 SECONDARY CMB ANISOTROPIES 

The solution of the Boltzmann equation for the present value 
of the perturbation of the photon temperature {6T/T)cmb, 
along the line of sight (los) 7, can be written as: 

(7)= to / ^x(x,r;)7- v(x,77), (1) 

where r) = /(I + z)dt is the conformal time and To = 
nefiO-TTioc, with Uefi bciug the present free electron density, 
ut the Thomson cross section and c the light speed. The 
quantities x(x, rj) and v(x, q) arc the ionization fraction and 
the peculiar velocity in units of c, respectively, calculated at 
position X = 7(770 — rj) and conformal time 77. 

In principle, a map of temperature anisotropics can be 
simply obtained by integrating eq. (1) along each los pass- 
ing trough random slices of the simulation boxes. However, 
the periodic simulation boundary conditions would artifi- 
cially enhance the anisotropy signal by a non-negligible fac- 
tor (Gncdin & Jaffc 2001). To prevent this spurious effect, 
we randomly flip and transpose each simulation box around 
any of its six edges, hence breaking the fictitious correlations 
introduced by the computational method. We consider 30 
(65) simulation outputs from z = 18.7 to complete reion- 
ization, i.e. z ~ 13 (2: ~ 8) for the L20 (S5) model. The 
output redshifts are optimized to completely cover the path 
between the initial and final redshift of the simulation. Al- 
though this method might now somewhat underestimate the 
true anisotropy signal as we miss the contribution of scales 
larger than the box, the results constitute a solid lower limit 
to such quantity. In addition, we emphasize that the size of 
our box (L = 20/i~^ Mpc) is one of the largest used up to 
now for reionization studies, hence making the large-scale 
missing power a less severe effect. 

The dimension of the map is set by the angle subtended 
by the simulation box at the highest redshift; for the adopted 
cosmology Omax ~ 9.25 arcmin. The spatial information on 
the ionization fraction is obtained from the radiative transfer 
sirrmlations, whereas the peculiar velocity field is provided 
by the A^-body sirrmlation. We repeat the above integral 
tion for a random realization with the same volume-averaged 
value of the ionization fraction x- A map of the temperature 
fluctuations due to ionized patches (i.e. the inhomogcncous 
part) is derived by subtracting the two maps. The result is 
shown in Fig. 2 for the S5 (top left panel) and the L20 (top 
right panel) model. 
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3.1 Anisotropy distribution 



The statistics of temperature anisotropy can be analyzed in 
terms of spherical harmonics, Yim- 
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The angular power spectrum, Ci, is then defined as: 
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There is a strict relation between the probed angular 
scale and the multipole / in the formula above, 9 ~ 180/Z 
degrees. Therefore the extension and resolution of our maps 
sets an interval in I in which our analysis is meaningful, i.e. 
4000 < I < 1.67 • 10'', corresponding to 1/3 of the map and 
the pixel scale, respectively. To analyze the maps and ob- 
tain the angular power spectrum we use the software pack- 
age HEALPixt(G6rski, Hivon & Wandelt 1999). The results 
arc shown in Fig. 1 together with the primary CMB power 
spectrum from WMAP data fitting (Spergel et al. 2003). 
The signal due to patchy reionization dominates the primary 
CMB power spectrum for I ^ 4000 and reaches a maximum 
amplitude of ^ 1.6 x 10~^^. The amplitude in the two mod- 
els is comparable. The power spectrum obtained here is in 
agreement with that derived by Gnedin & Jaffe (2001), and 
it is roughly an order of magnitude smaller than the one cal- 
culated by Santos et al. (2003) via a semi-analytical model. 
A last aspect which is worth commenting is the fact that the 
anisotropy keeps a rather flat level up to the highest signif- 
icant multipoles. That is the indication that the secondary 
anisotropy from reionization keep their structure at least up 
to the arcsecond scale. 



3.2 Comments on observability 

The detection of the signal from patchy reionization requires 
high sensitivity experiments that can reach large multipole 
numbers, since the peak of the power spectrum is expected 
to be at / of the order of few xlO^. These characteristics are 
within the capability of the next generation of millimeter 
wavelength interferometers like ALMA*, ACT§, or CQ^. 
For example, ALMA is expected to reach sensitivities of 
2 /xK rms for a 1' beam with 1 h integration up to 2 ar- 
cmin scale, thus appearing as a perfect instrument to search 
for signature of inhomogeneous reionization. 

However, to measure the power spectrum from patchy 
reionization, several other astrophysical signals must be 
cleaned out from the maps. In particular, the main fore- 
grounds in the angular range discussed here are the thermal 
Sunyacv-Zel'dovich (SZ) and the Poisson noise from faint 
point sources. Thermal SZ is expected to be negligible, at 
least after multifrequence cleaning, for observations at 217 
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Figure 1. Primary CMB power spectrum from WMAP data fit- 
ting (dotted line; Spergel et al. 2003), together with the angular 
power spectra from patchy reionization in the L20 (solid) and S5 
(dashed) model. 



GHz (Zhang, Pen & Trac 2004). More important is the fore- 
ground from unresolved IR and radio sources, which is sev- 
eral orders of magnitude above the reionization signal at 
I ^ few X lO''. Luckily, this foreground contamination can 
be described in terms of a simple power-law (White & Ma- 
jumdar 2003). Thus, a foreground measurement at I ~ 10* 
would allow to extrapolate its value at lower /, where it can 
be subtracted to obtain a clean reionization signal. For this 
reason, multifrequency observations are particularly suited 
to subtract such foreground contamination. For this tech- 
nique to be successful though, a good knowledge of the in- 
strumental noise is required. 



4 21 CM RADIATION FROM NEUTRAL IGM 

The 21cm hyperfine transition of neutral hydrogen in the 
IGM provides a powerful probe to study the era of cosmo- 
logical reionization. In this paper we use the results from 
the numerical simulation of CM, that we briefly describe in 
this Section. We refer to the above paper for further details. 

The emission of the 21cm line is governed by the spin 
temperature, Ts- In the presence of a CMB radiation with 
TcMB ~ 2.725(1 + z) K, Ts quickly reaches thermal equi- 
librium with TcMB, and a mechanism is required that de- 
couples the two temperatures. While the spin-exchange col- 
lisions between hydrogen atoms are too inefficient for typ- 
ical IGM densities, Lya pumping contributes significantly 
by mixing the hyperfine levels of neutral hydrogen in its 
ground state via intermediate transitions to the 2p state. If 
a Lya background ^ 9 x 10"^^ (1 + z) ergs cm~^ s~^ Hz"-' 
sr~^ is present at redshift z, Lya pumping will efficiently 
decouple Ts from Tcmb. CM find that the diffuse flux of 
Lya photons produced by the same sources responsible for 
the IGM reionization, satisfies the above requirement from 
z « 20 to the time of complete reionization. As the IGM 
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Figure 2. Top panels: map of CMB temperature fluctuations, (i5T/T)cmb i due to patchy reionization for the S5 (left panel) and the L20 
(right panel) model. The maps measure ^ 9.25 arcmin on a side. Lower panels: maps of differential antenna temperature, (5T/T)2icmi 
from 21cm line emission for the S5 (left panel) and the L20 (right panel) model, at an observed frequency (rcdshift) v 115 MHz 
(z 11) and V ^ 90 MHz (z 14.5), respectively. The maps have been derived assuming a bandwidth Av = 1 MHz and measure 
9.25 arcmin on a side. 



can be easily preheated by primordial sources of radiation 
(e.g. Madau, Meiksin & Rees 1997; Chen & Miralda-Escude 
2003), the universe will, most likely, be observable in 21cm 
emission at a level that is independent of the exact value 
of Tg. Variations in the density of neutral hydrogen (due to 
either inhomogeneities in the gas density or different ionized 
fraction) will appear as fluctuations of the sky brightness of 
this transition, and allow, in principle, to map the history 

of reionization 1 1 . 

Using the numerical simulations described in Section 2, 
CM have studied the evolution of the 21cm line emission ex- 
pected from those reionization histories. In particular, they 
have derived maps and fluctuations of brightness tempera- 
ture at different redshifts (i.e. observed frequencies) in both 
the S5 and L20 model. The S5 model predicts a peak in 
the amplitude of the expected rms brightness fluctuations 
at an observed frequency (redshift) v ~ 115 MHz {z « 11), 
whereas in the L20 the peak is at ~ 90 MHz {z « 14.5). In 
both models, the overall amplitude of the signal at its peak 
is {ST^^g)^^'^ « 10 — 20 mK at angular scales « 5 arcmin 
(see CM for details). In Fig. 2 maps of differential antenna 
temperature, {ST/T)2icm = 5Tt/TcMB{0), are shown for the 

I The brightness temperature against the CMB is defined as 
Tj, = TcMB^~^ + ~ e""^), where r is the optical depth 

of a patch of IGM in the hyperfine transition. The differen- 
tial antenna temperature between the patch and the CMB is 
Sn ^ (Ts -Tcmb)t/{1 + z). 



S5 (lower left panel) and the L20 (lower right panel) model 
at the peak frequencies. The maps have been derived from 
the simulations of reionization described in Section 2 assum- 
ing a bandwidth Av — uA.z/{l + z) = 1 MHz and measure 
~ 9.25 arcmin on a side. 



5 CMB/21 CM CROSS-CORRELATION 

CMB secondary anisotropies from patchy reionization are 
expected to be highly anti-correlated with 21cm line emis- 
sion temperature fluctuations on scales smaller than the an- 
gle subtended by typical Hii regions at the redshift of the 
21cm emission. To quantify this effect, we compute the cross- 
correlation between the CMB and the 21cm map at redshift 
z as: 

JVg L V -i / CMB \ I J 21cm J 

1=1 j=l 

where A^cmb (A^2icm) is the number of pixels in the CMB 
(21cm) map. The components of C{9,z) are then binned 
according to the separation angle 9 between the two lines of 
sight passing through the pixels i and j of the 21cm map at 
redshift z. No is the number of values falling in [6, 9-\-9min\, 
where 9min is the angular dimension of the cell. We repeat 
this procedure for a random binning. The results are shown 
in Fig. 3 and 4 for the S5 and L20 model, respectively. The 
subscript r refers to the random cross-correlation. The small 
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Figure 3. Cross-correlation between the CMB secondary 
anisotropy map and the 21cm maps at different redshifts (marked 
by different line styles) for a bandwidth Av = 1 MHz. In the 
small panel, a zoom of the region in which the transition from 
anti-correlation to correlation occurs is shown. The subscript r 
refers to the random cross-correlation. Model S5. 
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Figure 4. The same of Fig. 3, but for Model L20. Note that lines 
refer to different redshifts. 



box in the top-left corner shows a zoom of the region in 
which the majority of the curves first pass through a zero 
point, i.e. from anti-correlation to correlation. 

As expected, we find that the two signals are highly 
anti-correlated below a characteristic angular scale, 6o, ex- 
cept for the highest redshift where only a very small fraction 
of the volume is ionized (line at z — 18.5). The amplitude of 
the anti-correlation signal increases with decreasing redshift 



until reionization is almost completed. At any given redshift 
the model L20 shows a stronger anti-correlation signal, as 
reionization proceeds more rapidly, and a larger So, as a re- 
sult of the relative larger (on average) Hii region sizes. In 
fact, the angular scale 6o indicates the typical dimension of 
the Hii regions (including overlapping) at that redshift and 
allows, in principle, to reconstruct the reionization history 
and to discriminate among different reionization models and 
sources (e.g. quasars or massive Pop III stars versus more 
standard stars). 

The redshift evolution of the angular scale 8o is shown 
in Fig. 5; labels there indicate the typical comoving dimen- 
sion of the corresponding Hii regions in units of Mpc/h. As 
the evolution of 6o refiects the growth of Hii regions, the 
value of 00 at a given redshift is very different for the two 
models considered. In general, this result can be used to 
discriminate among different ionizing sources as Hii regions 
produced by quasars or massive Pop III stars typically tend 
to be larger than those digged by more standard stars. More- 
over, as the redshift evolution of 8o reflects the growth of 
Hii regions, non-monotonic reionization histories would re- 
sult in a more complex behavior for 9o. For example, in 
case of a double-reionization (e.g. Cen 2003; Wyithe & Loeb 
2003) we expect that 9o increases until the first reionization 
is completed. Then, once the ionizing emissivity drops and 
the IGM partially recombines, Oo should decrease or remain 
constant, and eventually grow again when the second reion- 
ization takes place. 

From our simulations it is also possible to derive a re- 
lation between the measured value of 9o and the volume 
averaged hydrogen ionization fraction in the computational 
volume. This relation is shown in Fig. 6. From there we see 
that the two quantities are positively correlated. It is also 
worth noticing that a 50% mean ionization level is char- 
acterized by 9o « 3' independently of the adopted reion- 
ization model: identifying this epoch is crucial as it corre- 
sponds to the redshift at which most of the CMB secondary 
anisotropics are produced and it provides a sensible defi- 
nition of the reionization redshift for prompt reionization 
models often adopted for practical purposes (Bruscoli, Fer- 
rara & Scannapieco 2002). Moreover, the Oq - Xe relation 
appears to be quite insensitive to the details of the reioniza- 
tion model, thus providing a robust mapping between the 
correlation function and the mean ionization level at each 
epoch. 

The results in Fig. 3 and 4 have been derived assuming 
a bandwidth for 21cm line observations of /S.u = 1 MHz. 
Although a smaller bandwidth, e.g. 0.1 MHz, would sub- 
stantially increase the intensity of 21cm line emission (CM), 
it does not significantly affect the estimates of the cross- 
correlation (Fig. 5). 

In conclusion, we find that the cross-correlation between 
secondary anisotropy in the CMB and 21cm emission maps 
can be a useful tool to follow the reionization process and to 
give constraints on the nature of the ionizing sources. More- 
over, the cross-correlation, combining information obtained 
by different experiments, can be used to maximize the signal 
from the reionization process with respect to instrumental 
noise, systematic errors in the measures, and astrophysical 
foregrounds. 
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Figure 5. Redshift evolution of the angular scale at which the 
transition between anti-correlation and correlation takes place in 

the L20 (upper curves) and S5 (lower) models. The labels report 
the typical comoving dimension of the corresponding Hll regions 
(including overlapping) in units of Mpc//i. The 21cm maps have 
been obtained assuming a bandwidth Ai/ = 1 MHz (solid lines) 
and 0.1 MHz (dotted). 
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Figure 6. Relation between the zero-point angle of the cross- 
correlation function and the mean ionization fraction in the 
simulation box for the S5 (left panel) and L20 (right) mod- 
els. The 21cm maps have been obtained assuming a bandwidth 
Ai/ = 1 MHz (solid hnes) and 0.1 MHz (dotted). 



6 SUMMARY 

We have calculated the secondary anisotropies in the CMB 

temperature power spectrum produced by inhomogeneous 
reionization from radiative transfer simulations consistent 
with WMAP observations. We find that the signal arising 
from this process dominates over the primary CMB com- 
ponent for / ^ 4000 and reaches a maximum amplitude of 
Z(Z -I- l)Ci/27r ~ 1.6 X 10"^'' on arcminutc angular scales, i.e. 
I as large as several thousands. We then cross-correlated the 
secondary CMB anisotropy maps with 21cm line emission 
fluctuation maps for the same reionization simulations. As 
expected, the two signals are highly anti-correlated on an- 
gular scales corresponding to the typical size of Hii regions 
(including overlapping) at the redshift of the 21cm map. 
The cross-correlation and, in particular, the redshift evolu- 
tion of the angular scale at which the transition between 
anti-correlation to correlation takes place, can be used: (a) 
to study the nature of the reionization sources, (b) to recon- 
struct the cosmic reionization history, (c) to infer the mean 
cosmic ionization level at any redshift. 

Cooray (2004) has studied the correlation signal be- 
tween CMB temperature anisotropies and 21cm fluctuations 
by means of an analytical reionization model. He concludes 
that, contrary to what we have shown, the correlation can- 
not be seen in the angular cross-power spectrum, due to a 
geometric cancellation effect between velocity and density 
fluctuations. However, such cancellation likely occurs due to 
the assumption made in that paper that the neutral hydro- 
gen fraction depends only on overdensity but not on spatial 
location. In fact, Cooray writes the neutral hydrogen density 
as uhi = XHn{l + 5), where xh is the neutral H fraction, 
n is the mean gas density and 5 is the gas overdensity. In 
a patchy reionization scenario, where the ionized bubbles 
around luminous sources do not completely fill the cosmic 
volume, it is clear that this assumption is not correct, as 
two fluctuations with the same 5 value located either inside 
a ionized region or outside it will have different uhi- Hence, 
if the patchiness of the reionization can be properly modelled 
{e.g. through radiative transfer simulations) the degeneracy 
(and the above cancellation) can be broken. Notice that the 
largest contribution to the secondary anisotropies come from 
the epoch when roughly 50% of the cosmic volume is filled 
with bubbles and where the variance in the relation nni — S 
is largest. 

Planned millimeter wavelength interferometers, like 
ALMA and ACT, are expected to have sensitivities and 
angular resolution good enough to measure the signature 
of inhomogeneous reionization in the CMB maps. However, 
extracting information on the reionization process from the 
observed maps can be hampered by the presence of both 
astrophysical foregrounds and instrumental noise. The same 
applies to 21cm emission observations (see Di Matteo, Ciardi 
& Miniati 2004 for a detailed study of the foreground con- 
tamination of 21cm maps). Provided that both the CMB 
and 21cm maps can be cleaned from foreground contami- 
nation, the information obtained from a cross-correlation of 
the two maps is an invaluable tool to study the reionization 
history and its sources. 
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